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Abstract 

Two suggested models of jet fragmentation in a quark-ghion plasma have been 
tested, combined and further developed. This has been done by generating hard pro- 
cesses in proton-proton collisions at -^/s = 2.76 GeV in PYTHIA 8 with a modified 
parton shower algorithm. Subsequently, a jet analysis of the final state hadrons have 
been performed with the FastJet recombination package. The results have been com- 
pared to unmodified jets in proton-proton collisions at the same collision energy, and 
to experimental data obtained by ATLAS. It has been shown that modifications of the 
parton splitting kernels alone is an insufficient modification to reproduce experimental 
data on dijet energy asymmetry and azimuthal decorrelation. Additional jet decollima- 
tion qualitatively reproduces the main features of the data provided certain variables 
are chosen aptly. We also study other jet characteristics with the model, such as the 
jet substructure, that may be observable experimentally in a near future. 



1 Introduction 



Quantum chromodynamics (QCD) is the theory of quarks, gluons (jointly referred to as 
partons) and their interactions. It bears similarities to quantum electrodynamics (QED), 
not least in the existence of charge. Whereas there is a single QED charge (electric), the 
QCD charge comes in three varieties referred to as colors (e.g. red, green and blue) and 
is carried by quarks. Similarly anti-quarks carry the corresponding anti-color. Gluons, of 
which there are eight, are the force particles of QCD. In contrast to photons, gluons are 
not charge neutral but carry a combination of color and anti-color, meaning that gluons 
are allowed to self- interact. 

Another remarkable difference between QED and QCD is the coupling, Og, which de- 
pends on the momentum scale used. At small momentum scales Og rapidly becomes large 
(confinement) whereas at large momenta it tends to zero (asymptotic freedom) [HIS]- The 
QED coupling has an opposite trend, however, can in most cases be considered constant. 

The concept of asymptotic freedom has led to the realization that at very high temper- 
atures or densities interacting partons become free and transform into a deconfined phase 
of matter - the quark-gluon plasma (QGP) OHj. The matter of the universe occupied this 
state a few microseconds after the Big Bang. At present not many details are known about 
the dynamics of this exotic state of matter. 

The appropriate conditions for the formation of a QGP can be obtained in a laboratory, 
namely in high energy collisions of heavy nuclei. This provides the possibility to study its 
properties in a controlled environment. Attempts in the 1980s and 1990s at SPS (CERN) 
led to the announcement of indirect evidence for the creation of the new state of matter 
[5]. Formation of a QGP was confirmed at RHIC (2005) [6j, where studies are currently 
carried out using Au-Au collisions at ^/s = 200 GeV per nucleon. The recent completion of 
the ALICE, ATLAS and CMS detectors in the Large Hadron Collider (CERN) has led to 
new data gathered from Pb-Pb collisions at ^/s = 2.76 GeV per nucleon, i.e. significantly 
higher than at RHIC. 

The data indicates that the bulk of the particle production originates from a hot source 
in thermal equilibrium. One of the most striking features is the so called "jet-quenching" 
phenomenon. This manifests itself in the suppression of particles with high transverse 
momenta at RHIC, as seen in Fig. [T| and the LHC. The yield of pions created in head-on 
heavy-ion collisions is suppressed by a factor ~ 5 compared to the expectation, had the 
nuclear collision simply been a mere superposition of proton-proton collisions. This implies 
that energetic partons transversing a dense QGP lose significant fractions of their energy, 
predominantly from radiative processes [71 [8] . 

The purpose of this investigation is to qualitatively reproduce some of the obtained 
data by constructing a basic model of the parton behavior in a QGP. Comparing it to 
similar data for proton-proton collisions, we thus hope to develop a better understanding 
for the plasma itself and the interactions between it and the particles traversing it. This 
discussion will mainly focus at the results obtained at the ATLAS experiment at the LHC. 
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Figure 1: Factor 5 supression of the vr'^ meson yield in heavy-ion colhsions. Data, taken 
from Ref. [9] , is compared to a model for radiative energy loss [30] . 



We will start with a brief discussion of the running of the strong coupling and how 
this results in the idea of a quark-gluon plasma. Then follows an outline of the underlying 
theory and mathematics relevant to factorization of hard processes, culminating in the 
DGLAP equation. We then proceed to discuss the importance of jets; jet fragmentation; 
and hadronization in the Lund model, and how to define jets for further implementation. 
The main investigation then follows with a set of reference results for intrajet distributions, 
energy distributions and jet energy asymmetry obtained for proton-proton collisions in a 
vacuum. We then introduce two ideas for modifications that are believed to reproduce 
data on jets in a QGP obtained at the LHC. The previous distributions are obtained again 
with the new modifications, followed by a discussion of parameter values relevant to the 
models. We finish the investigation with a summary of the results and what conclusions 
can be drawn from them. 



2 Theory 

2.1 The running coupling and perturbative QCD 

There are two main ways of solving QCD: on the lattice; and using the perturbative 
approach. The former involves discretizing time and space and considering the values 
of the quark and gluon fields at all the vertices of the resulting 4-dimensional lattice. 
This approach is suitable for calculating static quantities in QCD. However, the following 
discussion will consider high energy systems where lattice QCD is not suitable due to the 
large number of lattice points required for a complete solution. Instead perturbative QCD 
will be advised. 
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Figure 2: The strong coupling as runs with the momentum scale Q in an inverse logarithmic 
manner. Figure provided by Ref. [ID] 



Perturbative QCD is the idea of an orderly expansion in a small coupling = |j <C 1, 
where gs is the strong coupling constant. Some observable / can be calculated as 

f = fo + has + fial + hal + ... (1) 

where a fixed number of terms only need to be considered as the remaining ones should be 
negligible, assuming that the /j-terms do not compensate the smallness of a^. 

In high order QCD calculations it is often convenient to introduce a renormalization 
scale ;U, in order to preserve consistent dimensions for all quantities. The coupling constant 
Us depends on the scale at which it is evaluated: 

= , , ^2 , as{Q^) = , ^ (2) 
6o In J— bo In 

In eq. ^ the coupling is expressed in terms of the renormalization scale and in terms 
of a more conventional momentum scale Q applied to collider events. Here Aqcd is a 
constant scale at which the coupling diverges, and for which, if // 3> A (or as{fi'^) <C 1), 
perturbative QCD is valid, and 6o is (almost) a constant [Ijp]. It is suggested from eq. 
Q that at low momentum scales quarks are strongly bound, i.e. confined into hadrons, 
whereas at high momentum scales the coupling becomes weak and the quarks behave as if 
they were asymptotically free. 

As mentioned previously, asymptotic freedom has led to the idea that as some suffi- 
ciently high momentum scale the strong coupling should be weak enough to allow a decon- 
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Figure 3: Cross section for some arbitrary hard process qq — t- ab 



fined state of partons. It is sometimes useful to express the strong couphng as asiksT). 
Various calculations have been done in the attempt to find the critical transition temper- 
ature where the partons become deconfined jll^ [T2] . This is an example of a calculation 
where QCD on the lattice is a suitable approach, and is also employed by Ref. |12j. It 
is believed that the transition temperature is of the order 0(100) MeV corresponding to 
more than a trillion Kelvin. This temperature range would have been reached a few mi- 
croseconds after the Big Bang, after which the resulting quark-gluon plasma would have 
cooled down below the critical temperature and hadronization would have begun, forming 
the commonplace hadronic matter of every day life. 



2.2 Factorization properties 

In predicting the probability of any event to occur it is common practice to consider the 
cross section, a, of the event. In the simplest of cases one may consider two incoming 
particles in a head on collision resulting in two outgoing particles. The cross section of 
such an event would depend on the momentum fractions of the incoming particles, Xi 
{i = 1,2), carried by the colliding quarks, and the total energy transfer in the collision. 
For processes at high energies, it factorizes into the following form, cf. Fig. [3] namely 

CTft = y dxi/g/p(xi) j dX2fq/pix2)cr(^gg^ab)ixiX2s). (3) 

Here s is the square of the center of mass energy such that the partonic cross section 
depends on the energy transfer (3:1X25). The subscript h has been added to emphasize 
that the process is hard. The function fq/p{xi) is the number density of quarks of type q 
carrying a fraction xi of the momentum pi of one of the incoming baryons (and similarly 
fq/p{x2) for the other one). These functions are referred to as parton distribution functions 
(PDFs) and encode the long-distance dynamics of the process. The partonic cross section, 
cr(qq^ab) i ou the othcr hand describes the short-distance process. 

Quarks being charged with color may well end up radiating a gluon at some stage in 
the event. It is thus of importance to consider what corrections for cj/j would be required 
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Figure 4: (left) Initial state parton splitting, (right) Virtual term. 



in such an occasion. The following line of argument will consider only the case of an initial 
state emission (Fig. [4| left), however it should be noted that equivalent conclusions would 
be reached if final state emission was regarded instead. 

The correction to the cross section is certainly going to depend on the momentum 
fraction of the incoming quark, z, carried by the emitted gluon. However, an additional 
complication enters the picture in that the probability of emission depends on the momen- 
tum (resolution) scale, Q, of the process. 

. , , .UsCf dz dk'^ ^ iV^-1 4 
a,,.{p)^a.{zp) — — ^, C. = ^ = - (4) 

Here z is the momentum fraction the quark attains from the incoming parton, Nc = 3 
is the number of QCD colors, and kt = {1 — z)psm9 is the transverse momentum of the 
gluon with respect to the axis of the incoming beam. In the first term after the equating 
sign in eq. Q it should be noted that the hard cross section is now dependent on the 
remaining momentum of the quark, zp. The following term describes the strength of the 
interaction. The penultimate term relates to the momentum fraction, (1 — z) carried by 
the gluon. Note that the cross section diverges for z — )• 1. This is an important result and 
suggests that the gluon should carry a comparably small momentum i.e. the gluon is soft. 

The final term of eq. Q describes the momentum carried by the gluon in a direction 
perpendicular to the original axis. It is customary to employ transverse variables to avoid 
biases from boosts in colliders. Notice that the expression diverges for kt — )• 0, thus sug- 
gesting that the emitted gluon should carry very little transverse momentum i.e. emission 
should occur at a small angle compared to the beam axis. This is referred to as collinear 
emission. 

Certainly it is not wished for to have cross sections diverging to infinity. It turns out 
that if virtual terms are included, such as Fig. |4] (right), the soft divergence is cancelled 
out. To resolve the issue of collinear divergence a lower integration limit, ^p, is introduced. 
This cutoff, referred to as a factorization scale, creates a limit kt < hf below which all 
gluon emission are absorbed by the PDF from eq.Q. As a consequence the PDF becomes 
dependent on the cutoff (Fig. [5]) and the total correction to the cross section may be 
written 



5 




Figure 5: The correction to the hard cross section depends on the momentum scale Q and 
the factorization scale ^p- 



(Th+g+v - ^^^^ ^ / t^^*^^^^) ~ (^hixp)] q{x, fip) (5) 

where V represents the virtual term. The longitudinal momentum fraction x of the quark 
originating from the incoming baryon has been included. All gluon emission oikt < fip are 
absorbed by the PDF and exclusively those of larger transverse momentum are explicit in 
the 0{as) term in eq. (|5j). This term is now finite, however possibly large if Q ^ fip- 

2.3 DGLAP evolution 

Allowing the PDFs to absorb soft particles suggest the possibility of renormalization. The 
result is the DGLAP equation [HI El US]. The equation describes the evolution of the 
PDF as a function of the factorization scale: 

where Pqq is the real part of the 'Altarelli-Parisi splitting kernel' (i.e. a measure of the 
probability) for a quark to radiate a soft gluon. The + subscript signals that the virtual 
part is absorbed by the expression such that divergence cancels as z — ?• 1. Eq. ^ describes 
how, when the factorization scale is increased, extra partons with momentum fraction x 
are observed. These come from the branching of partons at lower factorization scales but 
with higher momentum fractions x/z (Fig. [6]). 

The incoming proton consists not only of quarks but of gluons as well. For the complete 
picture all possible evolutions must be coupled into the DGLAP equation: 
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Figure 6: DGLAP evolution, (left) State of quarks and anti-quarks only at = 12 GeV^ 
evolves into a state of quark-gluon mixture as increases to 150 GeV^ (middle left), 
(middle right) Purely gluonic state evolves into quark-gluon mixture at higher momentum 
scales (right). Figure provided by Ref. [T6] 



where the ® symbol is a shorthand notation for the integral in eq. ([6]). 

In this discussion only Pqq and Pgg will be considered, describing the emission of a 
soft gluon from a quark or another gluon respectively. It turns out that in the soft and 
collinear limits the two splitting kernels become equivalently dependent on z. It will thus 
be convenient to assume a common kernel for the two processes of the form 

It is worth noting that the coefficients Cj? = | and Ca = 3 reveals that gluons radiate 
other gluons with roughly twice the probability of quarks. 



2.4 Jets 

Jets are crucial to collider physics since they pose an important test of perturbative QCD. 
Normally, the fragmentation process of parton 6, not explicitly shown in Fig. |3| has 
to be considered. This is governed by long distance dynamics and is described by the 
DGLAP equations [T7j, analogously to the PDFs. Thus jets become a bridge between 
the measurable particles in the detector and the partonic subprocess. In particular, while 
the latter hard process takes place on short timescales and is usually assumed not to be 
modified by the presence of a deconfined medium, jets traversing a hot plasma are expected 
to be considerably affected. 



2.4.1 Jet fragmentaiton 

Returning to the simplest of situations in Fig. [s] (section 2.2) the event qq — )• ab may, for 
the purpose of argument, be interpreted as a. qq ^ qq type event. The (anti)quarks are 
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Figure 7: The jet must not be sensitive to soft and/or coUinear emissions, or hadronization 
effects. Figure provided by Ref. |T6] 




hard so gluon emission is expected. From eq. ([4]) the gluon is hkely to be soft and close 
in angle to the emitter. Certainly the gluon may split into further partons and so on, and 
if sufficient energy is available the situation may turn into a parton shower where each 
succession carries a smaller and smaller fraction of the available energy. As it seems one 
could end up with an indefinite spreading of particles from the initially so simple set up. 
However, it turns out that a parton emitted at some stage in a parton shower must always 
be of smaller angle than partons in a previous succession, thus resulting in increasingly soft 
and collinear emissions |18j . This process is referred to as angular ordering. Rather than 
considering these particles individually it is convenient to cluster them into cone shaped 
jets of some angular resolution. The jets should not be sensitive to the softest and most 
collinear emissions but remain unaffected by such happenings. The process is shown in 
Fig. [7} Notice that if a hard gluon is emitted it may well result in a jet of its own. 



2.4.2 Hadronization - the Lund model 

In nature, quarks are not seen as free states but are confined into hadrons. The process 
of hadronization is not trivial to calculate perturbatively as the hadron production must 
indeed involve a very large number of gluon interactions, which due to their softness is 
not suitable for a perturbative approach (c/ eq. ([T]) and eq. section 2.1). However a 



phenomenological model (the Lund model) reproducing the main features of hadronization 
has been put forward [19j . 

Similar to electromagnetic charges, QCD color charges have associated field lines be- 
tween all (anti)quarks and gluons. These lines tend to collapse into a 'string' of strong color 
field, with a fixed radius Ar ~ 0(1) fm, hence storing a roughly constant amount of energy 
per unit length. Thus, in the event of quark separation a long flux tube is produced. If a 
qq pair is produced within the tube and therefore removing a section, then the remaining 
energy stored within the tube is reduced. 

Due to the force along the fiux tube the hadrons will easily pick up longitudinal mo- 
mentum, however any transverse component remains small (0(200) MeV/ c). Thus jets 
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(a) Color field lines collapse 
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Figure 8: The process of jet formation according to the Lund model. Figure provided by 
Ref. m 




Figure 9: Jet production in colhder events |T6 



should be insensitive to most hadronization in addition to soft and collinear emission (Fig. 
[7| (right)). 



2.4.3 Jet definitions 

In some collider events it may be obvious how to define the jets from the outgoing particles 
such as the dijet event in Fig. [9] (left). However, considering the (identical) center and 
right events it is not trivial to deduce whether three or four jets are observed. In order to 
do so a precise way of defining a jet must be provided. 

A jet definition is essentially a set of rules by which a computer can take a list of 
particle momenta for some event, and cluster them into jets. These rules must obey the 
insensitivity to soft and collinear emission. The type of definitions considered here are 
known as recombination algorithms. They work in a backwards manner by intending to 
invert the splitting of the parton shower and recombining the particles. This does not only 
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provide a set of jets but also the possibility of iterating the process and inspect the contents 
of the jet. Note that this procedure also reconstructs so-called "fake jets", i.e. random 
clusters of particles. 

The three recombination algorithms considered here take the form 
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(9) 



The idea behind such a set up is to find a measure for the distance, dij, between every pair 
of particles This is then compared with some maximum allowance called the beam 
distance, diB- If the inter-particle distance is the smaller one the particles are 'combined' 
into a single new particle - a pseudojet. On the other hand, if the beam distance is 
the smaller one, then particle (or pseudojet) i is denoted a jet and removed from the 
list of particles. The process continues until no particles remain. In eq. ^ kt is the 
transverse momentum with respect to the collision axis (normally the z-axis). R is the 
opening (half) angle of a cone within which the momentum vectors of the particles in the 
jet are located. ARij is the distance between particles i and j in (j) — r] space such that 
Ai??j = {(pi — + (t]i — r/2)^, where 4> is the azimuthal angle of the particle and rj is 
the pseudorapidity defined as 77^ = — In [tan (g)] with 9 being the polar angle between the 
collision axis and the particle momentum. 

The three algorithms differ in the order they cluster the particles. The kt algorithm 
starts with the softest particles whereas the anti-/ct algorithm begins with the hardest. The 
latter ensures that the jets are centered around a hard 'core particle', evolving in concentric 
circles, thus forming jet cones of well defined shape. The Cambridge/ Aachen algorithm 
does not consider the transverse momentum at all but relies exclusively on the angular 
separation of the particles. 



3 Jets in a vacuum 



With the theory in place it is now time to start the implementation. The objective of 
this discussion, as previously stated, is to test, and further develop, models suggested to 
replicate data obtained for hard processes taking place inside a QCD medium. However, 
before any QGP investigations are carried out it is essential to analyze particle behavior 
in a vacuum. This will be done using proton-proton collisions at a center of mass (CM) 
energy of 2.76 TeV. This energy is chosen as it corresponds to the CM energy of heavy-ion 
collisions at the LHC performed in 2010. The investigations are carried out using event 
generator PYTHIA 8 |21j . which we have used to generate events with a hard transverse 
momentum scale pt = 70 GeV/c. Subsequently, the final state particles are reconstructed 
using the Fast Jet package (2.4.2) |22| . which implements the recombination algorithms 
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mentioned in section 2.4.3 A minimum cut of 20 GeV is placed on all measured jets in 
order to exclude fake jets. 

PYTHIA 8 is a full fledged event generator, incorporating as one of its key features a 
Monte Carlo parton shower algorithm. Using the Altareli-Parisi splitting kernel and the 
physical phase space, see eq. Q, it calculates the probability of a gluon emission with 
given z and kt and, thus, can generate a full cascade. The process continues until the 
ultimate gluons are below some cutoff scale where perturbation is no longer valid. This 
produces a parton shower. Realistic events surely consists of hadrons rather than partons, 
thus hadronization is considered by PYTHIA 8. This is implemented using the previously 
outlined Lund model. 

Fast Jet allows for an appropriate comparison between the three relevant recombination 
algorithms in order to select one for further use. A useful investigation to carry out is that 
of intrajet distributions, i.e. how does the number of particles, N, within the jet, with a 
specific value of a certain observable vary as the value of the observable changes? In the 
following, only the hardest jet in the event is considered, and the constituent particles are 
measured. The variables investigated are: 

• ^ = — log (j-^^^^^Y^^ , where ^ measures the fraction of jet momentum pj^t contained 
in the constituent particle. 



cos ^ 



Pparticle'Pjet 
PparticlePj et 



, gives the angle between the particle and the main jet, and 



Pt = PpartideS'mO, which gives the particle momentum component perpendicular to 
the main jet, here referred to as the transverse momentum (note the difference to kt 
which in this discussion is used in reference to the beam axis). 



The fragmentation functions are illustrated in Fig. 10 , A few things can be said about 
the jet as a whole, and the different jet definitions. Starting with ^ it is evident that the 
distribution forms a hump and is hence referred to as the hump-backed plateau. It presents 
a mean of ~ 4.3, corresponding to a momentum fraction of ~ 1.4%. It thus appears as if 
the jet is composed of a collection of soft particles rather than a few hard ones. This is a 
well known distribution that has been theoretically predicted from perturbative QCD and 
is a signal of angular ordering [231 [211 ES]- The three algorithms agree fairly well for ^ < 3.4 
where the C/A and anti-Zcj algorithms peak. This corresponds to a larger momentum frac- 
tion of ~ 4.0%. It appears that the kt algorithm has a large sensitivity to soft particles. 
This is no surprise as the algorithm uses a 'bottom up' type of approach, clustering the 
jet constituents starting from the softest particles to the hardest ones. This leads to jets 
without well defined shape - a feature which becomes even more evident when looking at 
the distribution over 9. Despite defining a 'cone' of angle R within which the particles 
should be, the kt algorithm ends up with larger angles than the allowed R = 0.4. The kt 
'cone' thus becomes larger and every jet contains more soft particles, which also explains 
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Figure 10: Intrajet distributions in vacuum for the kt (black), C/A (blue) and anti-A;t (red) 
recombination algorithms. R = 0.4. 
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the larger values for the three distributions. It is worth mentioning that the number of 
jet constituents with ^ < 2 and ^ > 6 corresponding to momentum fractions > 12% and 
< 0.3% respectively are very rare, underlining the fact that very soft and very hard jet 
constituents are not common. 

The ^-distribution demonstrates that the anti-fc^ algorithm alone provides jets of a reg- 
ular cone shape. This is clear from the 9-axis intercept at precisely 0.4 for anti-fcj but larger 
values for the other two algorithms. This is a useful trait as it allows precise geometrical 
limits to be put on the jets, hence this algorithm is selected for further use in this project. 
Non-regular jet boundaries has practical implications, which is often held against the kt 
and C/A algorithms. An example is in the case when an area of a detector is misbehaving. 
It is not trivial to see how far a kt jet must be from that area not to be influenced by it. 
This and other reasons have led to the anti-/ct algorithm being adopted as the default jet 
algorithm at both ATLAS and CMS. 

It appears that the jet constituents primarily have momenta in a direction close to that 
of the jet as seen by the peaks of the ^-distributions occurring at 0.05 — 0.10 radians. This 
is as expected from the general collimated structure and formation of jets as was discussed 
in section [2.4.11 

The p(-distribution indicates again that the majority of jet constituents carry a very 
small momentum compared to the jet as a whole. The distribution appears to obey an ex- 
ponential law up until values of ~ 2 GeV/c after which the number of recorded readings 
become very small and the statistics unreliable. 

A further insight into the distribution of energy within the jet can be done by defining 
a transverse momentum density p±{r) at some angle r < R relative to the axis of the jet. 
Then the fraction of transverse jet momentum relative to the beam axis, contained within 
the 'subcone' of opening angle r is 

, . Jo P±{r')dr' 

^(^) = — ?ii — (10) 

Pt 



Fig. 11 illustrates what fraction of the jet energy is contained within the subcone. The 
jets appear rather collimated with approximately 75% of the energy concentrated in the 
innermost half of the jet and the most central third carrying about 50% {R = 0.7). It 
is worth noting that the more energetic jets appear more collimated (blue) than the less 
energetic ones (red). This means that despite increased radiation, the harder jets still carry 
most energy in the central core particles of the jet. 

It is of interest to study how the jets resulting from a collision balance each other. A 
convenient way of doing this is to measure the jet energy asymmetry, Aj, and the azimuthal 



angle separation between the two hardest jets (Fig. 12), defined as 
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Figure 11: ~ 75% of the jet energy is concentrated in the innermost half of the jet. More 
energetic jets are more collimated than less energetic ones. 



Aj = t^\l^\ \^4>\ = \<i^i-H. (11) 

Here Eti and Et2 are the transverse energies of the hardest and second to hardest jets 
respectively, and |A(/)| is the angle between them. Thus, an asymmetry close to zero in- 
dicates that a collision has resulted in two jets of approximately equal energies. A value 
strongly deviating from zero indicates that more than two jets are present. 

Here a limit have been put on the jet transverse momentum in order to avoid contri- 
butions from background from the underlying event. For the hardest jet Eti > 100 GeV 
and for the second jet Et2 > 25 GeV. Additionally, we demand that both of the jets are at 
mid-rapidity (|7y| < 1). The asymmetry distribution in Fig. 12 suggests that a large cone 
results in many jets whereas a small one results in a few. This is largely a result of the 
cutoffs described above, introduced to reduce the contribution from soft background par- 
tons. Thus, some jets, clustered with a small value for the parameter R, are not centered 
around a hard core particle (such as may be the case in Fig. [o] (right), section 2.4.3) and 
will be eliminated. 
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Figure 12: Asymmetry and back-to-back angle distributions in vacuum for anti-/cj algo- 
rithm. = 20 GeV/c. 



It appears that if a large jet cone is used, then most events will result in two jets 
whereas for very small cones the number of jets per event increases, as is seen from the 
progressively shallower slopes at lower opening angles. This is certainly expected due to 
the sheer geometry of the problem. It is worth noting is that a very large cone may engulf 
several particle clusters into one jet, thus miss a lot of interesting physics. It will be con- 
venient for further reference to limit the cone size to 0.4 < R < 0.7. 

The I A(/) I -distribution indicates that most often the two hardest jets are separated by 
an angle close to vr radians, again supporting the idea that most events result in two major 
jets, possibly accompanied by one or a few soft ones. 



4 Jets in a QCD medium 

Experiments at RHIC have established a strong suppression of hadrons with high transverse 
momentum in hard heavy- ion processes [261 EZ] • This suggests that there is an interaction 
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Figure 13: In a QGP medium induced gluon emission should result in additional soft 
partons on the expense of hard hadrons. During the formation time r needed for the gluon 
to decohere from the initial parton shower a typical number of collisions N^-^h = \ occurs 
giving the gluon a transverse 'kick'. 




between high pt partons and QCD matter resulting in a significant loss of energy for the 
partons in addition to the vacuum fragmentation - a process referred to as 'jet-quenching'. 
It has been suggested that the underlying reason is medium-induced gluon radiation [28] . 

The introduction of a background medium requires the mathematical models of QCD to 
be modified to account for two main effects: Primarily all hard particles must be prompted 
to radiate gluons at a higher rate than in vacuum. Secondly, the increased gluon emission 
must cause a momentum spread among the particles. 

A suggested model is that the partonic scattering cross sections are dominated by 
small-angle scattering. Multiple small-angle scatterings has the consequence that all com- 
ponents of the parton shower undergo Brownian motion, a trait that can be used to quantify 
the spreading of the jet energy [29]. 



Fig. 13 pictures the emission of a gluon from some hard particle during the forma- 
tion time r. It is straight forward to calculate the momentum uncertainty, hence using 
Heisenberg's uncertainty principle (HUP) find an expression for r: 

Ap = ^'(l - z)2p2 +p2 _^ I ^2p2 _^ p2 Pt_ 

P (12) 

AtAp = tAp ~1 =^ r~^ = ^ 

Pt PI 

where w is the energy of the hard particle. Notice that in a vacuum pt = usinO ~ lvO 
meaning that t oc ^ for a fixed angle. Thus in a vacuum a soft particle takes longer to 
decohere from the parton shower than a hard one does. 



Returning to Fig. 13, during the gluon formation time, in a medium the particle may 
undergo a series of collisions. The typical number of collisions can be expressed as Ncoh = j 
where A is the mean free path of the particle. If on average each collision contributes 
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with transverse 'kick' ^ the average transverse momentum, {pi) = iJ^Ncoh = ir^- It is 

Pt 

convenient to denote by q the average squared transverse momentum per unit path length, 
that the gluon accumulates in the medium. Thus; 

2 n:- 1^ 



Pt= yqi^, y^' 

from which it is evident that in a medium the soft particles decohere faster than the hard 
ones, i.e. contrary to the behavior in a vacuum. Thus, in this model of a QGP softer 
gluons decohere faster from the parton shower due to small-angle scattering induced by 
the medium. 

The intensity of the gluon emission across the medium of size L can be expressed as 

,^ as f Cf\ dio L 
al oc 



/CF\duj qL ^ ' 



a. 



2iT \CaJ ^ Vq^ 
otg ( Cf\ duj [to, 
2Tr[cA " 



where in the final line some critical energy scale, ojc, has been identified due to the dimen- 
sionless nature of the intensity. The color factor Cp (Ca) applies for radiation off a quark 
(gluon). By inspection ooc = qlP', leading to the observation that 




t{uJc) = \ — = L, (wc) = qL (15) 



i.e. the formation time for the gluon is the entire size of the quark-gluon plasma, meaning 
that the critical energy previously identified in fact is the maximum energy the gluon is al- 
lowed to receive from the medium. Thus an upper limit is introduced to the momentum gain 
of the partons, hence also a limit on the momentum spread. It is clear that for soft gluons 
of energy co < \/qL will be completely decorrelated from the initial jet direction, thus the 
observation is consistent with the idea that soft particles decohere more rapidly than hard 
ones. 



4.1 Toy model for qualitative observations 

Based on the above considerations, we proceed with a simple implementation of two main 
medium modifications of the jet fragmentation that we hope will capture the main physical 
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mechanisms involved. Firstly, the spectrum in eq. (14) suggests an enhancement of soft 



gluon production due to interactions with the medium. Secondly, soft gluons accumulate 



transverse momentum along their path through the medium, see eq. ( 15 ), resulting in their 
Brownian motion. This, in turn, will lead to a decoUimation of the soft constituents of a 
given jet. 

Our approach to implementing the former of the above outlined effects is a modification 



of the Altarelli-Parisi splitting kernels (eq. ([8]), section 2.3), suggested by Borghini and 
Wiedemann [30j . The idea behind this modification is to significantly increase the likelihood 
for soft gluon emission by enhancing the soft part of the splitting function as follows 

pgq ^ f^^\ 2 + fmed 

\CaJ l-z ^ ^ 

where fmed is a constant. At this stage in the investigation it is merely of interest to 
observe what qualitative eff'ects the modification may have. Therefore, fmed is arbitrarily 
set to equal one. Borghini and Wiedemann [30] suggests that the modification should be 
constant, fmed = 0.8. This value is chosen as it has allowed the aforementioned to replicate 
the suppression of the vr" meson in Au-Au collisions at ^/~s = 200 GeV, see Fig. [l] 

An idea presented by Casalderrey-Solana, Milhano and Wiedemann [29j is to assign 
a 'kick' to all particles traversing the quark-gluon plasma, a concept referred to as "jet 
collimation" . We assume that the induced gluon emission also takes place purely perpen- 
dicular to the beam axis. Every final state hadron is assigned an additional momentum 



kick of the order vgly perpendicular to their trajectory, see eq. (15). In our implemen- 
tation, this extra boost follows a Gaussian distribution with standard deviation a = qL 
in accordance with the above argument. For the moment, L is is assumed to depend 
linearly on the azimuthal angle, (p. Again, qualitative eff'ects are sought for so q is ar- 
bitrarily set to 1 GeV^/fm. The path length through the medium is allowed to vary 
according to L = L* + (Lq — L*)^, where L* is the maximum allowed medium size at 
^ = 0, and Lo is some lower cut. Here L* = 1.0 fm and Lq = 0.5 fm. Thus it is 
assumed that the shortest distance through the medium occurs at zero angle in azimuth. 

The ATLAS and CMS Collaborations reported after the first ^/s = 2.76 TeV Pb-Pb 
collisions in 2010 their observations of jet energy asymmetry and azimuthal angle separation 



between the hardest and second jet ^[8j. The relevant results are shown in Fig. 14a As 
before, a limit has been put on the jet transverse momentum in order to avoid contributions 
from background from the underlying event, namely E^i > 100 GeV and Et2 > 25 GeV. 
The data indicates an increased level of asymmetry for the Pb-Pb collisions compared to 
proton-proton collisions in vacuum. A clear broadening is seen; a shift to a higher mean; 
and a shift of the zero peak to a larger value. Here a cone opening angle of R = 0.4 
was used. The ATLAS Collaboration further reported that the asymmetry increased for a 
smaller R. 



It is evident from Fig. 14b that merely modifying the splitting kernel with a constant 
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Jet energy asymmetry 




(a) ATLAS data (b) Modified splitting kernel (c) Additional jet collimation 



Figure 14: 

(a) Asymmetry and back-to-back angle of the two hardest jets as measured at ATLAS. 
Yellow bins are PYTHIA events for pp collisions in vacuum at ^/S = 7 TeV, empty circles 
are corresponding experimental data and solid black circles are data from Pb-Pb collisions 
at \/5 = 2.76 TeV. R = 0.4. 

(b) Results from modified Altarelli-Parisi splitting kernel. Solid line shows vacuum results 
and dotted line shows medium results. 

(c) Results for modified splitting kernel in combination with jet collimation. 

fmed f 

q = 1 GeVVfm 
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does not replicate the characteristics observed at ATLAS. The result indicates that this 
rather alters the asymmetry and azimuthal separation such that they behave as if a smaller 
opening angle R was used. Fewer entries in the histogram suggests that the jets have less 
energy. This supports the proposition that gluon emission has increased and somehow 
transports energy out from the jet cone. However, the general behavior of the jets remain. 
This is increasingly obvious if one considers Fig. [Tsj It is clear that the number of particles 
involved in the hard process increases as the splitting kernel is modified. The ^-distribution 
reveals that the number of very hard particles, i.e. those of low ^, have decreased somewhat 
whereas the softer particles have increased in numbers. This indicates that the hard core 
particles of the jet are suppressed by the medium and radiating its energy as soft gluons as 
expected. The shift of the peak in the ^-distribution suggest a slight spread of the particles, 
however any alteration in the pt distribution is insignificant. It may thus be concluded that 
the main effect of this modification is to increase to soft emissions. 



Considering Fig. 14c , picturing the results from the collimation modification, it is clear 
that the behavior is dependent on the jet cone opening angle, R. It must be noted that the 
particles in the jets are assigned with additional transverse energy of GeV/c. For a 
large jet cone (blue) this additional momentum may not be enough to remove most gluons 
from the jet but rather spread them to the edges of the cone, thus resulting in an overall 
increased broadened transverse jet energy distribution. However, for a small jet cone (red) 
the decoherence may be quite sufficient to strip the jet of most gluons, hence reducing the 
over all energy. This explains the significant increase (reduction) of entries in the histogram 
for R = 0.6 {R = 0.2). For R = 0.4 (green) these effects seem to approximately cancel 
and the number of entries are of the same order for the vacuum (solid) and the medium 
(dotted) models. Just as the ATLAS and CMS Collaborations observed [71|8], in Fig. 14 a 
clear broadening is seen in the asymmetry; a shift to a higher mean, here by approximately 
30%; and a shift of the zero peak to a larger value which seems to increase somewhat 
for smaller jet cones. The A(/)-distribution reveals a tailing towards low values starting at 
~ 2.5 for R = 0.4, in rough agreement with ATLAS and CMS. 

It appears that a random kick in the transverse plane [29] in combination with the mod- 
ified splitting kernel |30j provides the correct qualitative features of the jet energy asym- 
metry distribution and the back-to-back behavior of the hardest and second jet. However, 
for a more realistic picture a discussion of the geometry of the collision, and appropriate 
selections of the parameters fmed and q are in place. 

4.2 Refined model - Nuclear geometry 

At the energies dealt with the incoming lead nuclei will be Lorenz contracted and appear 
virtually two dimensional in the collision moment. In a simple model the nuclei can thus 
be assumed to take the shape of circular 'pancakes' of radius Rnuc — 



with A = 207 being the mass number for lead. Fig. 16a illustrates how the overlap of these 



pancakes forms an almond shape in a semi-peripheral collision with impact parameter 
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Figure 15: Intrajet distributions for processes in vacuum (red) and in quark-gluon plasma 
(blue) {R = 0.4). 



b = 3.2j^ For the time being it will be assumed that this shape is filled with a quark-gluon 
plasma from the time of the collision, and remains static and unchanged throughout all 
measurements. It is necessary to calculate the distance from a jet creation point p(xo, Vo) 
to the edge of the almond. Assuming the path through the medium is a straight line and 
that the jet exits through a segment of a circle centered at {xao, VAo) one can write; 

XA = XAo + Rnuc COS 9 X = Xo + ^ COS (j) 

yA = VAo + Rnuc sine y = yo + ishKp (17) 

for circle A for linear path 

^The impact parameter is defined as the distance between the center of the two coliiding nuclei. 
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(a) Geometry of Pb-Pb collision (b) Random creation points and path length 



Figure 16: 

(a) Simple geometrical representation of a collision between lead nuclei, labelled A and B. 

(b) (top) Randomly generated points of creation for jets, (bottom) Distribution of path 
lengths through quark-gluon plasma. 



22 



and similar for circle B. Equating gives 



^ = [Rnuc ~ i^y COS (j) — sin (f))'^^^ ^^"^ — Ax cos ^ — Ay sin ^ (18) 

where Ax and Ay are the coordinates of the jet production point relative to the center of 
the pancake through which edge the jet exits. ^ is now the distance to the almond edge as 
a function of the jet direction in the transverse plane, the azimuthal angle and the point 
of production. Thus a randomly selected production point can be assigned to all created 
jets as is illustrated in fig. 16b (top). Here the distance between the nuclear centers is 
3.2 fm corresponding to an overlap of roughly 70%. The bottom picture indicates that the 
majority of particles travels less than 8 fm through the medium with a mean around 5 fm. 

Using ^ as the path length L it is time to return to the previous model. There is no 
longer any need to assume a linear behavior of L in (p. L varies as derived above and it is 
reasonable to expect that the modification of the splitting kernel could vary linearly in L 
such that 

fmed{L) = Lf^^^ + /^g^ (19) 

where f^^^ is some minimum modification (corresponding to a jet production point at 
the edge of the almond overlap and traveling outwards, hence set to zero) and f^^^ is a 
characteristic modification. Since we have no handle on the space-time structure of the jet 
fragmentation, we will additionally assume that all particles originate from a point where 
the hard process has taken place. This calls for an improvement to be addressed in the 
future. 

We have tested two options: (i) a constant value /^^^ = 0.8 |30] and /^^^ = 0, and 
(ii) a variable fmed{L) with parameters f^^^ = and /^g^ = 0.8/5, due to the mean path 



length being close to 5 fm, see Fig. 16 



Fig. 17a immediately reveals that the difference between a variable and a constant 
fmed is insignificant. Here the histograms have been normalized to the number of entries 
to simplify the comparison with the ATLAS data (Fig. 14), and R = 0.4. It is clear that 
modifying the splitting kernel alone is an insufficient model to describe the changes a jet 
experiences as it passes through the quark-gluon plasma. Adding the jet collimation mech- 
anism to the model replicates the ATLAS data with improved accuracy as is shown in Fig. 



17b, There is a clear broadening in the asymmetry distribution and significant increase of 



hard jets with azimuthal angle separation strongly deviating from vr. The model does not 
replicate the clear peak seen in the asymmetry spectrum obtained at ATLAS. Neither does 
it reproduce the distinct tail observed in the A(;/)-spectrum at ATLAS but rather alters the 
gradient of the plot. An interesting feature is that for asymmetries above ~ 0.22 there is 
no distinct difference between the vacuum and the medium situations. This may be a crit- 
ical value where three jet events become increasingly frequent, as is suggested by the little 
hump in the vacuum distribution. There is a vague indication in the ATLAS data that the 
gradients at large asymmetries should be of the same order. Thus three jet events appear 
to be less affected by the medium than dijet events. This can of course be experimentally 
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Jet energy asymmetry 



Jet energy asymmetry 




(a) Modified splitting kernel alone (b) Additional jet coUimation 

Figure 17: Asymmetry and azimuthal angle separation for refined QCD model. 



investigated. 

With a model roughly replicating the main characteristics of the measured ATLAS 
data it is of interest to move on to unexplored territory and investigate the intrajet struc- 
ture resulting from the model. Starting off with the energy distribution within the jet, Fig. 
18 shows the expected consequences of gluon emission in vacuum and medium. Angular 
ordering ensures that increased gluon emission results in a more collimated jet. Thus, 
harder jets tend to be more collimated as was seen in Fig. 11 (section^. However, in a 



medium the gluon spreading transfers energy to the edges of the jet cone, or even outside 
the jet. Hence, jets in a QGP are expected to be less collimated than jets in a vacuum as 
is shown in the simplest of cases in Fig. 18c Note that no adjustment has been made to 



the angular ordering as implemented by PYTHIA 8. As hard jets are more likely to emit 
gluons it is expected that these are affected more than less energetic ones. The results are 
shown in Fig. [T9| 

The above analysis agrees well with the observations for the intermediate jet cone 
(R = 0.4). For the most narrow cone {R = 0.2) it is clear that the medium has made the 
jets less collimated. However, it appears that the reduction in ip{r) is roughly the same for 
all energy ranges. This is probably related to the fact that most jets are almost completely 
stripped off of their soft constituents resulting in merely a small spreading of energy inside 
the cone. The largest jet cone {R = 0.7) presents more bewildering results. It is clear that 
the harder jets are feel the medium effects more extensively than softer ones (however not 
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(a) Angular ordering ensures addi- (b) The QGP spreads gluons towards (c) CoUimation for a jets in ( 
tional gluon emissions coUimates the the cone edge or even outside the jet. (black) and (b) (red), 
jet in vacuum. 

Figure 18: Analysis of jets consisting of a single hard particle emitting a gluon suggests 
that jets in a QGP should be less collimated than jets in a vacuum. 




(a) R = 0.2 (b) R = 0.4 (c) R = 0.7 

Figure 19: Jet collimation for various jet cones. Vacuum (solid), Medium (dotted). 
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much more than in the case of i? = 0.4). However, it appears as if the softer jets are less 
affected by the medium than they are for smaller cones. This may be an intrinsic artifact 
of the model, or possibly related to soft emissions from one cone entering the a neighboring 
one. 

Certainly the medium model produces rather large changes compared to the vacuum 
situation. As an example; considering the hardest jets in the central picture, the central 
half of the cone carries ~ 75% of the jet energy in vacuum whilst only ~ 35% in the 
medium. This poses the question whether the magnitude of the transverse momentum 
boost is of a reasonable size, i.e. what is a realistic value for the parameter ql 







— rr \ green:q = 1.0 
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Figure 20: Intrajet distributions for various values of q. 



The intrajet distributions in Fig. 20 demonstrate the differences between various values 
for q. The number of particles in the jet is strongly dependent on the transverse momen- 
tum assigned to them. For q = 0.16 GeV^/fm (purple) the number of jet constituents have 
increased compared to the vacuum situation (red). Thus, the modified splitting kernel ap- 
pears to dominate. However, for q = 1.0 GeV^/fm the number of particles have decreased, 
thus suggesting that most soft gluons have been pushed outside the jet cone. The modified 
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splitting kernel shows its effect in the reduction of particles with very low ^ (hight pt), 
whereas the jet collimation demonstrates its presence by reducing the number of particles 
with very high ^ (soft particles). The precise changes in these features would have to be 
determined experimentally and would reveal how the two effects balance. 

The ^-distribution reveals significant differences in the angular distribution of particles 
within the jet. For q = 1.0 GeV^/fm and q = 0.5 GeV^/fm the distribution forms a rather 
smeared out hump centered around ~ 0.16 and ~ 0.19 respectively. For q = 0.16 GeV^/fm 
the distribution is almost flat with the exception of the extreme edges. This is likely due 
to the fact that more collisions may occur before the gluon is completely decorrelated from 
the jet. The jet thus take the shape of only a few hard particles at the center, surrounded 
by an almost even distribution of soft particles filling up the jet cone to the edge. The pt- 
distribution demonstrates the same effect. Though most momentum is still carried by the 
central core particle, this fraction is significantly reduced on behalf of the softer particles 
that are spread to larger angles. 

5 Conclusion 

This investigation has been aimed at studying how the effects of jet-quenching can be 
modeled qualitatively. It has been performed by scrutinizing the mathematical theories 
behind perturbative QCD in the standard model and testing these as reference results. 
A basic model involving the modification of the Altarelli-Parisi splitting kernels and jet 
collimation by transversely 'kicking' particles has been attempted and used as a frame for a 
more serious model accounting for the nuclear geometry. Finally a discussion of the effects 
of the parameters fmed ^-nd q has been included. 

The result is a set of distributions for jet energy asymmetry and azimuthal angular 
separation between the two hardest jets in an event has been obtained. The distributions 
roughly replicate data collected at the ATLAS experiment in 2010, however with a few 
discrepancies. Overlooking the defects, the model has then been applied to intrajet distri- 
butions and the energy distribution within the jet. The result is a jet which in a medium 
strongly suffers suppression of its hardest constituents on the behalf of soft gluon emis- 
sion. These gluons rapidly leave the jet cone leading to a balance between two effects: an 
increase in gluon production and a successive removal of the same gluons. The balance 
between these two competing effects is a subject for further theoretical and experimental 
investigations. Overall, the jet suffers a significant spread of its constituent particles over 
the entire jet cone. The foremost part of the momentum is centered closely around the 
jet axis, however the jet energy is increasingly smeared out over the jet cone resulting in 
a less collimated jet. Each event appears to result in jets of angular separations strongly 
deviating from the tt radians commonly observed in a vacuum. Dijet events experience 
a significant change in asymmetry compared to the vacuum situation. However, events 
involving more jets appear to be less affected by the medium. Thus we can conclude that 
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the introduction of a quark-gluon plasma not only alters the jet behavior but also distorts 
the intrinsic structure of the jets. 

The model presented in the investigation relies on a set of assumptions, primarily that 
the physics derived for vacuum situations apply in the medium as well. In particular, 
the fragmentation of the hard process is assumed to be similar to a vacuum. There is no 
theoretical predictions behind the values fmed and q, and no consideration of longitudinal 
momentum is included, which would lead to an inelastic energy loss and stopping of the 
jet fragments in the plasma. The model is based upon a classical treatment of the nuclear 
geometry and assumes that there is no cross section associated to the area outside the 
nuclear overlap. A more appropriate approach would be to convolve the Saxon- Woods po- 
tentials from the two nuclei. Further the model assumes a stationary, uniform quark-gluon 
plasma, through which all jets traverse linearly. A more realistic picture would involve an 
expansion of the plasma with the consequence of reduced temperature and density. This 
would require a discussion of the thermodynamic properties of the plasma, which has not 
been a part of this investigation. 

Despite its discrepancies, the model reproduces some important qualitative features of 
the ATLAS data. It is useful starting point in the understanding of the underlying reasons 
for jet properties in a QCD medium. Further, it predicts some interesting results outside 
of the present experimental scope, which would be equally important to confirm as to 
disprove. Thus the investigation may be considered successful. 

References 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 

[7] 
[8] 
[9] 
[10] 



D.J. Gross, F. Wilczek, Phys.Rev.Lett. 30 1343-1346 (1973) 

H.J. Politzer, Phys.Rev.Lett. 30 1346-1349 (1973) 

N. Cabibbo, G. Parisi, Phys.Lett.B 59 67-69 (1975) 

J. Collins, M.J. Perry, Phys.Rev.Lett. 34 1353-1356 (1975) 



http://newstate-matter.web.cern.ch/newstate-matter/Experiments.html 



I. Arsene, IG Bearden, D. Beavis, C. Besliu, B. Budick, H. Bggild, C. Chasman, CH 
Christensen, P. Christiansen, J. Cibor, et al. Nucl. Phys. A, 757(1-2): 127 (2005) 



G. Aad et al. [ATLAS Collaboration], |arXiv: 1011. 618 2]y2 [hep-ex] (2010) 

[CMS Collaboration] , lar Xiv:1102.1957| y2 [nucl-ex] (2011) 

B. A. Kniehl, G. Kramer, B. Potter, Nucl. Phys. B 582 514 (2000) 



S. Bethke, Eur.Phys.J. C64 689 (2009) |arXiv:0908.1135| [hep-ph]] 



28 



[11 

[12 
[13 

[14 
[15 
[16 

[17; 
[is; 

[19 

[2o; 

[21 

[22 
[23; 
[24 
[25 
[26 
[27 
[28 
[29 



R. Hagedorn, Nuovo Cimento Suppl. 3 147-186 (1965) 

F. Karsch, E. Laermann, A. Peikert, Nucl.Phys.B 605 579-599 (2001) 
V.N. Gribov, L, N. Lipatov. Sov.J.Nucl.Phys. 15 438 (1972) 

Yu. L. Dokshitzer. Sov.Phys. JETP, 46 641 (1977) 

G. Altarelli and G. Parisi. Nucl. Phys. B 126 298 (1977) 

Graphics by G. Salam. ,arXiv:1101.2878v 2 [hep-ph] (10 Jan 2011) 

R.K. Ellis, W.J Sterling, B.R. Webber, QCD and Collider Physics, ch. 5, Cambride 
Univertisy Press (1996) 

Yu.L. Dokshitzer, V.A. Khoze, A.H. Mueller, S.I. Troyan, Basics of Perturbative QCD, 
ch. 4, Editions Frontiers (1991) 

B. Andersson, G. Gustafson, G. Ingelman, T. Sjostrand, Phys. Rep. 97 (1983) 

Graphics by G. Davies. Imperial College, Lectures on Nucl. and Particle physics (2011) 

T. Sjostrand, S. Mrenna and P. Skands, Comput. Phys. Comm. 178 (2008) 85 
|arXiv:0710.3820| 



M. Cacciari, G.P. Salam, G. SoyezPhys. Lett. B 641 (2006) | ;hep-ph/0512210] 

Yu. L. Dokshitzer, V. S. Fadin, V. A. Khoze, Phys. Lett. B 115 242 (1982) 

Yu. L. Dokshitzer, V. S. Fadin, V. A. Khoze, Z. Phys. C 15 325 (1982) 

A. Bassetto, M. Ciafaloni, G. Marchesini, A. H. Mueller, Nucl. Phys. B 207 189 (1982) 



K. Adcox et al. [PHENIX Cohaboration] , arXiv:nucl-ex/0410003 (2005) 



J. Adams et al. [STAR Cohaboration], |arXiv:nucl-ex/0501009| (2005) 

R. Baier, D. Schiff, B.C. Zakharov, Ann. Rev. Nucl. Part. Sci. 50, 37 (2000) 

J. Casalderrey-Solana, J.G. Milhano, U.A. Wiedemann, |arXiv: 1012.0745| i^l [hep-ph] 
(2010) 



[30] N. Borghini, U.A. Wiedemann, arXiv:hep-ph/0506218vl (2005) 



29 



